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Quantum Eletrodynamics (QED) is considered as the most successful of all physical theories. It 
can predict numerical values of physical quantities to a spectacular degree of accuracy. However, 
from the very early days it has been known that, in QED, there are two important problems which 
are linked with the very foundation of the theory. In 1952, Dyson put forward strong argumnts 
to suggest that the perturbation seires in quantum electrodynamics can not be convergent. Just 
three years latter, in 1955, Landau argued that the effective running coupling constant in QED has 
a pole (Landau singularity) albeit at some very high energy scale. This paper addresses, in details, 
the question of stability of perturbative vacuum state of QED in the light of these two well known 
problems. 

Landau has been a cult-like figure for many of us who studied theoretical physics in the former 
Soviet Union. As an undergraduate student in the department of theoretical physics of People's 
Friendship University, Moscow, in 1970's, I grew up hearing colourful stories about the legendary 
persona of Lev Davidovich from my teachers some of whom had known him at personal level. It is a 
great honour for me to contribute this article to the Landau centenary volume of Eletronics Journal 
of Theoretical Physics. 



Landau singularity of the effective coupling constant in quantum electrodynamics has been known for more than 
half a century [l|. The original derivation of the singularity was based on the summation of one loop diagrams of 
vacuum polarization tensor for photons in the perturbation theory. In the early days the validity of such an expansion 
was looked upon sceptically by many people including the authours themselves. However, after the advent of l/N 
expansion technique by t' Hooft, it was soon realised that this singularity appears at the leading order in 1/iV/ 
expansion where Nf is the number of specieses of electrons, also called the number of flavours. This implies that in 
the infinite flavour limit this singularity is exact provided the perturbation series in l/iVy expansion converges. This 
issue of convergence of the perturbation series will be one of the central themes latter in the paper. There have been 
attempts to interpret this singularity in many ways. Landau and Pomeranchuk tried to argue that this singularity 
reflects the fact that at short distances strong vacuum polarization effects screen the electric charge completely 
Others, including Shirkov, have called it Landau ghost reflecting the internal inconsistency of quantum electrodynamics 

It should be mentioned here that the stability of ground state and the possible existence of an ultra-violet fixed 
point has been studied extensively in the lattice formulation of QED in the wide range of value of the fermion flavour 
Nf by Kogut et al [1, 0|- Using lattice formulation, Kim et al has argued for the triviality of QED On the other 
hand in massless QED in the continuum, Miransky argues that there exists a chiral symmetry breaking phase [||. 
However, in these studies Landau singularity plays no role. There has been some recent studies of Landau singularity 
using lattice formulation of QED by Gockeler et al (Q and references there in). These studies seem to suggest that 
chiral symmetry breaking allows QED to escape Landau singularity. But then chiral symmetry breaking, as their 
study shows, seems to be intimately connected with trivility of QED. Landau singularity has also been considered 
from a different perspective by Gies and Jaeckel and by Langfeld et al @. The details of these approaches can 
be found in the cited references and will not be discussed here. Our approach will be very different from the ones 
described in the publications above. We will rather be interested in finding the meaning of Landau singularity than 
finding a way to escape it. 

We want to consider the issue of stability of "ground state/ vacuum state" of quantum theory in the unified and broad 
perspective of quantum field theory and the many-body theory. Therefore, at first, we consider many-body ground 
state of two purely quantum mechanical systems: "Coulomb system with large number of flavours of fermions" and 
"System of weakly interacting electron gas in a condensed matter system". These are discussed in section-I and 
section-II. After this, we return to the main theme of the paper. In section-Ill, we reproduce Dyson's argument for 
the divergence of perturbation series in QED. In section- IV, we show the connection between the Landau singularity 
and instability of vacuum state in QED. It is argued in the Dysonian framework, how the divergence of the series 
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removes both these problem. In section-V, we explain how a divergent asymptotic series can give rise to physically 
meaningful physical quantities. In the last and concluding section, we comment on the non-perturbative aspects of 
QED. Dyson's original arguments, as well as, our studies of Landau singularity and vacuum state shows that the 
physical observables in non-perturbative QED are non-analytic in the coupling , e 2 , as well as the inverse flavour, 
1/Nf, and the perturbative power series in these parameters can not capture this behaviour. In the absence of non- 
perturbative theory, we suggest that there should be attempts to experimentally search for non-analyticity in some 
physical observables in QED. 



I. GROUND STATE OF COULOMB SYSTEM WITH LARGE NUMBER OF FLAVOURS OF 

FERMIONS 



We will be investigating in this paper the question of stability of vacuum state in quantum electrodynamics. 
In this context, it is interesting to have some information regarding the stability of quantum machanical ground 
state of a large system of charged particle. Many-body theory of Coulomb system of fermions have been studied 
extensively in several publications. In quantum theory a many-body system with ground state energy Eq(N) is called 
thermodynamically stable or simply stable if Eq(N)/N is bounded below when the number of particles, N — > oo . 
In this context, the question of stability of matter consisting of negatively charged electrons and positively charged 
nuclei is very important. It was in 1967 when Dyson and Lenard proved that, in the framework of nonrelativitic 
quantum mechanics, matter consisting of N electrons and K static point nuclei of charge Z (= N/K) is stable 
10]. Subsequently, Lieb and his collaborators have made a very detailed investigation of the stability of matter in 
nonrelativistic as well as relativistic case ( [lH [l2| and references there in ). These studies seem to suggest that 
at high enough energies quantum eletrodynamics may not be a well behaved theory. Thermodynamic stability of a 
system of particles interacting via Coulomb interaction is associated with the control of the short distance behaviour 
of the interaction. In the nonrelativistic limit, zero point kinetic energy of the fermions controls this short distance 
behavoiur. However, in the relativistic limit, there is a need for certain bounds on the value of the fine structure 
constant [Hi [12) [13] • Landau singularity is also believed to be associated with the "high energy / short distance" 
behaviour of the electromagnetic interaction 

Let us consider the large flavour case in some detail. In 3-dimensional space, let us consider a volume of linear 
dimension R where there are N number of positively charged and N number of negatively charged fermions. 
Fermions of both positive and negative charges come with Nf flavours. Particles are assumed to be of the same mass. 
We assume that N >> Nf , however, both the number of particles and number of flavours are assumed to be large. 
Taking into account the the Pauli exclusion principle for the case involoving multiflavour fermions in 3-dimensional 

space, the kinetic energy K, apart from some numerical factors, can be written [l4j as K w lm j^n 2 / 3 ' potential 



energy U due to Coulomb interaction is given by U = — clq 6 ^ / , ao > is a numerical constant. The total 

energy, E(R), is given by 

The ground state is obtained by minimizing the total energy, E(R), with respect to R. 

ma 2 ., o I/a,., , ma?,, n i/^o 

E « ~—^(e 2 Ny 3 ) 2 N ; e = E a /N ^ - -/(e 2 iV. 1/3 ) 2 
h n 

h 2 N 1/3 1 R h 2 1 

R ° ~ ^a~ { W ^iVp ; r ° = {N/N f y/z ~ ma 0e 2iV / 1/3 (2) 

When the number of flavour Nf increases, the size of the small box tq decreases and when it approaches Compton 
wave length Ti/mc , we can no longer use the nonrelativistic quantum mechanics. Therefore, the above estimates are 
correct only when 

Nf « ~~7T ^ (3) 

In the relativistic quantum theory of many-body systems, it is almost customary now to take c\p\ fllj ]. where c is 
the velocity of light, as the kinetic energy of individual particles. Therefore, the kinetic energy K for the system of 
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N particles with Nf flavour in a region with linear dimension R is given by, K — (hcN A ^ 3 ) / (RNy 3 ) . Thus, in 
the relativistic case, 

E(R) = K + U « ^Ljj -a £ -^— (4) 

In this case stability of ground state requires E(R) > O.This leads to the condition, 

N f < ~^ — \ — (5) 

which in the unit c =Ti = 1 , becomes Nf < ( ao \syi ■ At this ponit a few comments are necessary. We considered 
a theory with large N and large Nf , and it is implicit in our discussion that we are considering, N — > oo , but we 
have not said much about the flavour Nf . From the equations above, we find that it is also possible to take the limit, 

Nf —> oo but slower than N. Assuming that in the limit Nf — > oo and e — > , e 2 Nl = constant, and that it 
satisfies equation Eq.(5), it is easy to conclude that, in this limit, matter is stable. This mathematical limit, has the 
physical implication that when the charge, e is small and the flavour Nf is large, the relevant parameter that sets the 

stability criteria is e 2 JVy^ 3 . It is the value of e 2 iV^ 3 that decides the stability of the system. When e 2 JV^ 3 > 4j , 
in other words, when the number of flavours Nf > , the many-body system is thermodynamically unstable. 

In subsequent section, we will find that in quantum electrodynamics the relevant parameter is not e 2 iV^ 3 but e 2 Nf. 



II. FELDMAN MODEL OF WEAKLY INTERACTING ELECTRON GAS 



The main theme of this paper, as annouced in the abstract, is to look for the meaning of Landau singularity of the 
effective coupling constant in QED. There are some other contexts in which the effective coupling constant develops 
singularity. We have in mind some condensed matter systems, where these singularities have well defined physical 
meaning. In this section, we refer to the renormalisation group analysis of weakly interacting Fermi system with short 
range potential at finite density and zero temparature by Feldman et al flU Ha . \vi\ . The iterative renormalization 
group transformations show that if there is an attractive interaction among the electrons in any angular momentum 
channel, then there appears similar type of singularity in some suitably defined running coupling constant. We briefly 
describe here the Feldman model of weakly interacting electron gas. This section is essentially a short description of 
the results taken from the very detailed review paper by Froehlich, Chen and Seifert [16(. 

The model of weakly interacting electron gas studied by Feldman et al is a condensed matter Fermi system in 
thermal equillibrium at some temparature T (for simplicity, assume T = 0) and chemical potential fi. On microscopic 
scale(w 10~ 8 cm), it can be described approximately in terms of non-relativistic electrons with short range two body 
interactions. The thermodynamic quantities such as conductivity depend only on physical properties of the system at 
mesoscopic length scale (rs 10~ 4 cm), and therefore, are determined from processes involving momenta of the order 
of zfi- around the Fermi surface, where the parameter, A >> 1, should be thought of as a ratio of meso-to-microscopic 
length scale. This is generally refered to as the scaling limit (large A, low frequencies) of the system. The most 
important observation of Feldman et. al. is that in the scaling limit, systems of non-relativistic (free) electrons in d 
spatial dimensions behave like a system of multi-flavoured relativistic chiral Dirac fermions in 1 + 1 dimensions. The 
number of flavours N w const. A . It is possible then to set up a renormalization group improved perturbation 
theory in 4 around the non-interacting electron gas, where in, the large number of flavours N, play an important role 
in actual calculations. 



A. Free Electron gas and the Multiflavour Relativistic Fermions in 1 + I Dimensions 

Let us consider a system of non-relativistic free electrons in d spatial dimensions with the Euclidean action, 

S {r,i>) =J2j d d+1 xrAx)(ido - - M )^(x) (6) 

The Euclidean free fermion Green's function, G° aa , (x—y), where a and a' are the spin indices, x = (t, x) and y = (s, y), 
t and s are imaginary times, t > s, is given by, 
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G% a , (x-y) = (C {xW„ (y)) „ = -<W f (dk) ^ lko(t ^ lk{x y) (7) 

Where we have used (dfc) = ( 2 ^)(d+i) d d+1 k. In the scaling limit, the leading contributions to G° aa ,(x — y) come 

from modes whose momenta are contained in a shell Sp^ of thickness around the Fermi surface Sp. In order to 
approximate the Green's function, let us introduce the new variables tD, py, po such that /cpu; G Sp, po = k and 

k = (kp +p\\)uj. If k € S^, then py << fcp, and we can approximate the integrand of Eq.(7), by dropping p| term 
in the denominator. In other words, 

G° aa ,(x -y) = <W J ^Lkp'e^^GS - s,uj(x-y) (8) 

where da(u) is the uniform measure on unit sphere and 

/-rfp rfp|| e -»fco(t-»)+i P ||J(i?-g) 

G c (t-s,w(x-y )) = - / — — ± : (9) 

J 2ir 2ir ip - v F p\\ 



is the Green's function of chiral Dirac fermion in 1 + 1 dimension, vp = kp jm is the Fermi velocity. The aJ-integration 
in Eq.(8) can be further approximated by replacing it with summation over discrete directions iJj by dividing the shell 
Sp^ into N small boxes B^.,j — l,..,N of roughly cubical shape. The box, B$., is centered at luj € Sp and has an 

approximate side length ^.The number of boxes, N = Qd-i^^ 1 , where fid-i is the surface volume of unit sphere 
in d spatial dimensions. The Green's function is, now, given by 



G° aa ,(x-y) = -S^J2 J 



dpo dp\\ p ± e -iPo(ts)+ip(x-y) 

2ir 2ir 2ir ipo — VfP\i 



(10) 



where p = pyu; + p± is a vector in — kpujj and po G 1Z. Thus in the scaling limit, the behaviour of a d- 
dimensional non-relastivistic free electron gas is described by N = ttd-i A d_1 flavours of free chiral Dirac fermions in 
1 + 1 dimensional space-time. The propagator G c (t — s,co(x — y) depends on the flavour index Q .But the energy of 
an electron or hole with momenta k , depends only on py, where py = k.u) — kp and u = .It is proportional to 
Py , just as for relativistic fermions in 1 + 1 dimensions. 



B. Renormalization Group Flow and the BCS Instability 



Large scale behaviour of the weakly interacting system is described by an effective action. To see how the effective 
action is calculated, let us consider a system with Euclidean action of the form, 

S(V*,V) = + W,^) (ii) 

where So(ip*,ip) is the quadratic part given in Eq.(3) and Si(ip*,ip) is the quartic interaction term. For a weakly 
interacting electron gas, 

S/(V*> VO = 9ok F - d J2 I d d+1 xd d+1 y : ^{x)^ a {x)v{x - $)S(x - yoWAvWAv) (12) 

cr,a> 

The factor kp~ d ensures that go is dimcnsionless. The two body potential v(x) is assumed to be smooth and short 
range, and therefore, its Fourier transform v(k) is also smooth. In order to calculate the effective actions, we first 
split the field variable into slow and fast modes, tp = -0< + ^> , where Supp ?/>> C R x (R d \ S^) (region >) and 
Supp C R x Sp' ' (region <) . We then integrate out the fast modes using the functional integral for fermions, 

e -Se// = I I e - s W>d<,i>>,i><) (13) 
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Now using the linked cluster theorem, we obtain 

e -s. ff (^<) = exp(-5 0: <- < Si > G o > +1 < Si; Si > G% ~ < Si; Si; Si > G o -...) (14) 

The abbreviations, < a; b > < a; 6; c > etc., denote the connected correlators. The subscripts G> indicate that 
the expectations < (.) > G o are calculated using infrared cut-off free propagators in accordance with the functional 
measure, 

dP = (l/E)D^* > D^ > e~ So ^>^>'> 

The connected correlators can be evaluated using the Feynman diagram technique. Therefore, the effective action 
can be calculated once we know the amplitudes of connected Feynman diagrams. It is clear that the S e ff contains 
far more interactions than the original quartic interaction Si. However, for weakly interacting systems the original 
coupling remains dominant. To carry out the iterative renormalization group scheme of Feldman et al, we choose some 
large scale Ao << ^= and calculate the effective action, S e ff perturbatively to leading order in -^-.The effective action 

depends on modes corresponding to wave vectors located in the shell, Sp , of width ^ around the Fermi surface. 

Now, we divide the shell Sp^ into N = const. A d_1 cubical boxes, B$ j , of approximate side length ^.Next, we rescale 

all the momenta so that, instead of belonging to the boxes Bq. they are contained in boxes B$. of side length w Uf- 
These two steps are generally known as decimation of degrees of freedom and rescaling. The renormalization group 
scheme consists of iteration of these two steps. 

Assume that the degrees of freedom corresponding to momenta not lying in Sp have been integrated out. Let 
^cr(fc), k G R x S^p \ denote these modes. The sector fields are defined as 

1>aA*) = I (dkJe^-^^Mk) (15) 

It is easy to see that ip a ( x ) = X)<3 elkFU3 ^^Ps.a{ x ) ■ Inserting the Fourier transform of sector fields in Eq.(6) and 
Eq.(12), and carrying out some algebraic manipulations, we obtain 

So = -J2f (^)^. ff b)(iPo-v F cJp + 0(^))^(p) (16) 

Si = Y k F~ d E v(kF(Zi -«3 4 )) / \dp 1 )...(dp 4 )(2ir) d + 1 S(p 1 +p 2 -p 3 - P4 ) 

Wi,.. ,uj^a,<7' 

^h u a(Pi)^ 2 ,a'(P2)^cs 3 M'{P3)^c3 i ,a{P4 : ) + terms of higher order in-j- (17) 

Using cluster expansions to integrate out the degrees of freedom corresponding to momenta outside the shell Sp®\ one 
can show that , at scale ^f- , the effective action has the form given by Eq.(13) and Eq.(14) except that v(kp{uji — w 4 )) 
is replaced by a coupling function g(uj\, (J2, W3, C04) ~ v(kp(uji — W4)) with uj\ + LU2 = ^3 + ^4 • 

Next, one considers the rescaling of the fields and the action. The fields are rescaled in such a that the supports of 
the Fourier transformed "sector fields" are boxes, Bj — X(B$ — fcfw), of roughly cubical shape with sides of length 
kF, and such that the quadratic part of the action remains unchanged to leading order in j. The first condition 
implies that p 1 — ► p = pX and x 1 — ► £ = § ■ The rescaled sector fields and their Fourier transforms are given by, 

foAO = A a Vto, CT (AO ; fe i<r (p) = X^-^Aj) (18) 

Inserting the scaled Fourier transformed fields into quadratic part of action So, it is easy to see that Sq remains 
unchanged if the scaling dimension a = |. Now, inserting the rescaled fields in the quartic part of the action, we find 
that the quartic part has scaling dimension (1 — d). The quartic terms of higher degree in momenta as well as terms 
of higher degree in fields appearing in the effective action have smaller scaling dimensions. Thus the effective action 
in terms of scaled sector fields is, 

(#)^, CT (p)(iPo-v F wp)^ i(T (p)+ S g(wi,w 2 ,W3,w 4 ) 

UJ,<J U?i+U72 = L03-\-LJ4-,(T.(T f 
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(dp 1 )....(dp 4 )(2TT) + 5{pi +p 2 -p 3 -P4)^S 1 ,aiPl)^a 2 ,a'iP2) 1 >Pu; 3 ,a'iP3)^uj i ,a(Pi) 

+terms of higher order in— (19) 

A 

We see that that the inverse propagator for the sector field is diagonal in 65 ,and it depends only on p and p|| = uop 
but not on p± = p — (65.p)65. 

We are interested in the renormalization group flow equations in the leading order in j. Therefore, for carrying out 
the decimation of degrees of freedom, we will be interested only in those diagrams that contribute to the amplitude 
in the leading order in j. Before we consider these diagrams, let us consider the possible inter sector scattering 
geometries. How many independent g^°\u5i,652, 65%, 65^) exists? For d = 3, suppose cJ 3 ^ — W4.0n the unit sphere, 
there are = Const. Aq -1 different 65's. But all choises of 65 1 and 65 2 with 65i + 65 2 = cD 3 + 65^ lie on a cone 

containing 65% and 65i with 65% + 65^ as the symmetry axis. Therefore, there are 0(Xq~ 2 ) choices. Only when 65% = —65^ 
that there are iV( ) = Const. Aq -1 choices. Couplings involving incoming states with uj 3 ^ —0J4 will be represented 
by g( Q \wi, U2, 65s, 604,). Couplings that involve sectors u5 3 = — or cquivalcntly 65 1 = — 65 2 will be denoted by 

ffscs^ 1 '^ 4 )- Because of rotational invariance, g^ cs (65\,6jA) is a function of only the angle between 65\ and 6J4. 

The chemical potential receives corrections from connected diagrams with two external electron lines (self energy 
correction of the electrons). The contribution of order zero in j- comes only from the tadpole and turtle diagrams. 
These diagrams contain one internal interaction squiggle of order }_ i and there are iV^ ' = Const. An -1 choices of 

the inner particle sector. Therefore, the contribution is of order zero in i-. The amplitude corresponding to these 
tadpole and turtle diagrams turns out to be p-independent but of order O(g^0)). Thus, Sfii = 0(g^ / Xq) , and the 
renormalized electron propagator is Gj 1 (po,p.65) = — [ipo — v-pp.u) + Ao<5/ii] _1 . 

To find the evolution of the coupling constant g(65i, 65 2 , 65%, W4), we have to calculate amplitude of diagrams with 
four external legs. It is found that when 65\ + 65 2 = ^3 + 6)4 7^ 0, the coupling functions do not flow in the leading 
order in j. But for sector indices 65\ + 65 2 = 6)3 + £4 = 0, the coupling functions, g^ cs (65\, 6J4), flow. The diagrams 
that contribute to the flow equation are the ladder diagrams with self energy insertion for the internal electron lines 
but with no other two legged subdiagram.The amplitude of such a diagram with n interaction squiggles and with zero 
incoming and outgoing box momenta of the particles is given by, 

(ti n , 65 n - 1 )...g B cs(65i, w') 

cJi ,...,uT n 

In the equation above, (65', —65') and (65, —65) are sector indices of incoming and outgoing electron lines respec- 
tively.Other sector indices correspond to the internal electron lines. [3 is a strictly positive number coming from the 
fermion loop integration in the Feynman diagram and is given by (3 = / dkj_dk\\dko [k$ + (vpk\\ — XS^ii) 2 }^ 1 . We 
find that the renormalized value of gscs is 

9^(65,65') = ^(65,65') + f^(-^) n £ (-l^^^ s (<3,<3„)^ s (<3 n ,a5„-l)...fl^ s ( ^ 3l,a5 , ) 



71=1 

7 0) 



+0( 9 —) (20) 

The explicit expression for flow equation can be obtained by expanding the coupling functions gBCs(65,65') = 
9bcs((65,6j') into spherical harmonics, gBcs((65,65') — ^ gihi({65,65'). Up to terms of order j, 

9 — i£f + <> (21) 

To obtain the differential equation for the R.G. flow, let us define g\^ := g[ X3 \ and consider a scale A = e'A . Next, 

define gi (t):=g\ etXo) . The coefficient (3 = /3(t, t') vanishes in the limit t' \ t, therefore, j3(t' , t) = (f - t)j(t) + 0((t' - 

t) 2 ). Writing difference equation for the couplings g^ +1 ^ and gj*\ and dividing both sides of the difference equation 
by (t' — t) , we finally obtain 

j t gi(t) = - 1 g l (t)(t) 2 + 0(e- t g(t) 2 ) (22) 
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where 7 = j(t) > . It is independent of I and approximately independent of t , and therefore, we set 7 = 70- The 
positivity of 7 follows from slow monotone growth of j3(i! — t) in t' . Neglecting the error term, the solution can be 
written as, 

#W = t~ tttt;- 23 ) 

If the coupling constants are positive or rather non-negative, gi(0) > 0, the effective running coupling constant goes 
to zero, and we have the Landau-Fermi liquid phase. This phase consists of free quasi-particles which are electrons 
and holes with renormalized mass and the chemical potential. On the other hand, if there is an angular momentum 
channel, I , with attractive interactions (gi(Q) < 0) the flow diverges at a finite value of the scaling parameter, 
t = — (7o<?;(0)) _1 ■ This singularity reflects the instability of the ground state. 

The ground state of the perturbation theory described above was taken to be the non-interacting Fermi gas with Fock 
space constructed from the elementary excitations, electrons and holes. The renormalization group analysis shows 
that this is not the true ground state in the presence of atrractive interaction. The singularity in the running coupling 
constant reflects just this fact. To see whether the true ground state is a BCS state, we need to know the nature 
of the pole. It is easily found that the residue at the pole of the effective running coupling constant has negative 
sign. This signifies the presence of Cooper pairs [1, [l8j]. The true ground state is thus, the BCS ground state of 
superconductivity. In the following section, we shall see that the Landau pole, which looks very similar to the pole of 
the effective coupling constant in condensed matter system, has positive residue at the pole. This sets it apart from 
the BCS type of instability @,[lj]. 



III. DYSON'S ARGUMENTS FOR THE DIVERGENCE OF PERTURBATION SERIES IN COUPLING 

CONSTANT 

Dyson's arguments for the divergence of perturbation theory in QED is elegant in its' simplicity. We will simply 
reproduce here his arguments [l9j |. 
Let us suppose that 

F(e 2 ) = a + aie 2 + a 2 e 4 + ... 

is a physical quantity which is calculated as a formal power series in e 2 by integrating the equations of motion of 
the theory over a finite or infinite time. Suppose that the series converges for some positive small value of e 2 ; this 
implies that F(e 2 ) is an analytic function of e at e = 0. Then for sufficiently small value of e, F(— e 2 ) will also 
be a well-behaved analytic function with a convergent power series expansion. However, for F(— e 2 ) we can also 
make a physical interpretation. In the fictitious world, like charges attract each other. The potential between static 
charges, in the classical limit of large distances and large number of elementary charges, will be just the Coulomb 
potential with the sign reversed. But it is clear that in the fictitious world the vacuum state as ordinarily defined is 
not the state of lowest energy. By creating a large number N of electron-positron pairs, bringing the electrons in one 
region of space and the positrons in another separate region, it is easy to construct a pathological state in which the 
negative potential energy of the Coulomb forces is much greater than the total rest energy and the kinetic energy of 
the particles. Suppose that in the fictitious world the state of the system is known at a certain time to be an ordinary 
physical state with only a few particles present. There is a high potential barrier separating the physical state from 
the pathological state of equal energy. However, because of the quantum mechanical tunneling effect, there will always 
be a finite probability that in any finite time-interval the system will find itself in a pathological state. Thus every 
physical state is unstable against the spontaneous creation of many particles. Further, a system once in a pathological 
state will not remain steady; there will be rapid creation of more and more particles, an explosive disintegration of the 
vacuum by spontaneous polarization. In these circumstances it is impossible that the integratation of the equation of 
motion of the theory over any finite or infinite time interval, starting from a given state of the fictitious world, should 
lead to well-defined analytic functions. Therefore F(— e 2 ) can not be analytic and the series can not be convergent. 
The central idea in Dyson's arguments for the divergence of perturbation theory in coupling constant, as is evidient 
from the lengthy discussion above, is that the convergence of the perturbation theory in coupling constant would 
lead to the existence of pathological states to which the normal states of QED would decay. These pathological states 
correspond to states of a quantum field theory whose vacuum state is unstable. Therefore, if quantum electrodynamics 
is a meaningful theory, the perturbation series must diverge (for more discussions related to these arguments, see 
[20L l2l| and references there in ) . 

It should be noted that Dyson's proof appeared much before the advent of asymptotic freedom in quantum field 
theories. The main point in Dyson's proof is that, if the perturbative series is convergent, then for small value of 
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e 2 , we can analytically continue to — e 2 and then this series will also be convergent. Let us consider the series for 
the vacuum polarization (two point Green's function for photons). Both the perturbative as well as the analytically 
continued series are assumed to be converegent. We can carry out loop wise summation. We can write the formal 
sum for one loop diagrams ( the coupling is assumed to be small) and extract from it the effective coupling constant. 
In the analytically continued theory, we obtain 

e 2 

^ = 1 - (-e 2 )ln£ (24) 

When A — > oo, then e 2 ^ — > , and therefore the analytically continued theory is asymptotically free. It is also easy 
to infer that at low energies the effective coupling constant increases. On purely formal grounds, the efective coupling 
constant of this purely formal theory behaves in the same way as the effective coupling constant in QCD for both the 
high energy and low energy limits inspite of the fact that in the formal theory interactions are mediated by abelian 
gauge fields. For further considerations, we require an infrared regulator. We can always choose an infrared regulator 
such that the effective coupling constant remains small. It is not hard to see now that the classical argument of 
Dyson holds. Let us consider a box of size L with wave functions vanishing at the boundary. The size, L, itself can 
be taken as the infrared regulator. One may suspect that, as the the number of pairs of electron-positron created, as 
per Dyson's arguments, increases, the inter-particle distance decreases and this leads to decrease of e 2 ^, and since 
e 2 ^ goes to zero as the inter-particle distance goes to zero ( because of asymptotic freedom) , the vacuum state 
somehow stabilizes. However, it turns out that the suspision is unfounded, and the reason behind this is that the 
effective coupling constant, e 2 ^, decreases too slowly with the decrease in the inter-particle distance. Let us suppose 
that TV-pairs of electron-positrons are created, the electrons separate to one half of the box and the positrons to the 
other half. In this process, the vacuum energy decreases at least by Eq ~ —N 2 e 2 ^j ( other factors are suppressed 
). Now e 2 e ff, for large N, goes as, e 2 ^ ~ l//n(7V 1 / 3 ). Therefore, the decrease in the vacuum energy by the process 

of pair creation goes as Eq ~ —N 2 e 2 e ^ = — N 2 /In^N 1 ' 3 ), which approaches — oo as N — + oo. This demonstrates 
that the vacuum energy remains unbounded below. Therefore, asymptotic freedom does not change the pathological 
character of the analytically continued theory and Dyson's argument remains intact. 



IV. LANDAU SINGULARITY AND THE VACUUM STATE 

The Lagrangian of QED with number of flavours, Nf, is given by 22], 

£ = J2^ 3 +m- erfAp) ^ + -F 2 „ (25) 

where ipi and rjj^ are the Dirac field and its' conjugate, j is the flavour index, and and F M „ are the electromagnetic 
potential and the field strength respectively. We will investigate this model when the number flavours, Nf, has both 
the positive and negative sign. We, therefore, introduce the notation \Nf\ = sign(Nf) x Nf . For the Lagrangian, 
written above, ^- expansion is introduced by assuming that, in the limit |iV/|— >oo, e 2 \Nf\ — constant = a 2 (say) 
. Large flavour limit of quantum electrodynamics can be equivalently described by the Lagrangian, 



£ = E ¥ (i-fd* +m- YA,y + -Fl v (26) 

With this form of the Lagrangian, it is easier to set up Feynman diagram technique. To each photon and fermion line 
corresponds their usual propagator. Each vertex contributes a factor of — y= = , each fermion loop contributes a factor 

of (—1) for anticommuting fermions and a factor of Nf because of summation over fermion flavours. Using these 
rules, it is easy to set up 1/Nf expansion series for any physical observable. Just as in the case of perturbation 
theory in the coupling constant, the expansion in allows us to express an observable F in the form, 

F ^ = «° + Wf Ql + N] Q2 + (2?) 

Qo , Qi , Q2, are some functions of the coupling contant. Now suppose that the series converges for some 

small value of j^j ( large value of Nf ), then the observable function F(-^j) is analytic for = ( Nf = 00 
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). Therefore, we can consider a small negative value of ( large negative value of Nf ) for which the function is 
analytic and convergent. In other words, the function F(jj-) can be analytically continued to small negative value 
of 1/Nf and the series thus obtained will be convergent. Latter, in the text, we will discuss the meaning of negative 
flavour. Here we just mention that, in the context of lattice QCD, fermions with finite number of negative flavours 
have been considered before(see [23,l2J,l23] and references there in). 

Let us calculate the effective coupling constant from the formal 1/Nf expansion series of the two point Green's 
function. The series is assumed to be convergent, and therefore, for sufficiently large Nf , one can restrict to the 
leading order term. The leading order term is given by the one-loop diagrams which can easily be eavaluated to obtain 
the polarization from which one can read off the effective coupling constant. It is given by, 



If Nf is positive, 



=e//(A 2 ) = f Nf (28) 

37r|iV/| m 2 



37r m 2 



It is to see that the effective coupling constant has a pole at finite but very large value of A 2 = m 2 exp(iir 2 / e 2 ) . 
This is known as the Landau singularity of the effective coupling constant in QED. This is the central theme of the 
paper, and therefore, we repeat that the convergence of the 1/Nf expansion series allows us to choose a sufficiently 
large Nf and restrict to the leading order terms in 1/Nf. 
Thus :- 

• The appearance of Landau singularity in the effective running coupling constant in QED is intimately linked to 
the assumption that the 1/Nf expansion series converges. 

We will now argue that Landau singularity leads to the instability of vacuum state in QED. We mentioned before 
that, the convergence of the series for sufficiently large positive Nf, allows us to analytically continue it to negative 
Nf, and the series will again be convergent. Thus, for Nf negative, we obtain, 

«5»< A '> = rr?bf (30) 

From this equation, it is easy to see that in the limit A — ► oo , e 2 ^- — > 0. Therefore, the formal theory that we obtain 
from the analytical continuation of 1/Nf ( for large Nf ) to the small negative value of 1/Nf , is ( at least formally 
) asymptotically free. This seem to suggest that the physical meaning of the negative sign of Nf could possibly be 
traced in the free theory without the interaction term. In following sections, we shall argue that the choice of negative 
Nf for anticommuting fermions amounts to considering commuting fermions with positive Nf. 

For the Lagrangian given by Eq.(25) (in four dimensional Euclidean space), the partition function is given by functional 
integral, 

Z ac = J DA(x)Di>(x)Dip(x)exp(- J d i x£) (31) 
Funtional integration with respect to the anticommuting fermion fields (grassman variables) gives, 

Z ac = J DA{x) det N f (iffy + m - erfAJ 

exp{-\ ( d A xF 2 ) (32) 



On the other hand, if we consider the fermion fields to be commuting variables, then the integration above ammounts 
to functional integration over complex fields, and we obtain, 

Z c = J DA(x) det~ N f (i^df, +m- e^A^) 

exp{~ J d^xF 2 ,) (33) 
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Note that this expression could be obtained from the previous expression, simply by assuming that Nf is nega- 
tive. Therefore, anticommuting fermions with negative Nf has the same partition function as the commuting fcrmions 
with positive Nf. Since, physically interesting observables can be calculated from the partition function, our claim is 
that the negative flavour anticommuting fcrmions are equivalent to the positive flavour commuting fermions. 

We can also arrive at this concision using the Feyman diagram tecnique of the formal perturbation theory. Consider 
the two point Green's functions for the photons using Lagrangian given by Eq.(26) . First, we consider just one loop 
diagram and show how the contribution due to flavours appears in the calculations. There are two vertices and a 
fermion loop, each vertex contributes a factor of = , the fermion loop contributes a multiplicative factor of 

(— 1) because the fermions anticommute and a multiplicative factor of Nf because of summation over flavours of the 
internal fermion lines. Now, if Nf happens to be negative, the factor (—1) and the factor Nf , combines to give the 
factor \Nf\ .This is also the contribution if the fermions commute and the flavour is positive ( the factor (—1) is 
absent for commuting fermions). The same procedure applies for the multiloop diagrams. This shows that the choice 
of negative Nf for anticommuting fermions amounts to considering commuting fermions with positive Nf . We 
have shown earlier in text that quantum electrodynamics with anticommuting fermions and negative value of Nf 
is asymptotically free. Therefore, our purely formal quantum electrodynamics with commuting fermions and positive 
Nf is asymptotically free [22l . It is well known that the quantum field theory of free commuting fermions does not 
have stable vacuum state [26l l27l [28j . It then follows that an interacting asymptotically free quantum field theory 
built around such a vacuum state can not be stable: all states in this theory will be pathological. These results follow 
from the single assumption that the l/7V/-expansion series in QED is convergent, and therefore analytic in 1/Nf. As 
per Dyson's argument, this would lead to the decay of normal states of QED with anticommuting fermions to the 
pathological states of QED with commuting fermions via the process of quantum mechanical tunnelling. Therefore, 
the vacuum state of QED with large number of flavour of anticommuting fermions can not be stable. 
We have already seen that the convergence of 1 /Nf expansion theory invariably leads to the appearance of Landau 
singularity. 
Thus :- 

• Landau singularity signals the instability of vacuum state of quantum electrodynamics. 

We have, no where, in text shown that the 1/Nf expansion series diverges. It was only an assumption. There 
exist large number of publications which, based on the behaviour of the large order terms in the series (in coupling 
constant), suggest that the perturbation series is, possibly, divergent. Similar arguments can be extended to our case 
unless there is some magical cancellations in large orders of the series. But such magical cancellations, if any, would 
plague the theory with the instability of vacuum state and render it meaningless. 
Thus :- 

• Quantum electrodynamics with large number of flavours of (anticommuting) fermions will be a meaningful 
theory, only if the series in 1/Nf expansion diverges. 

In the abstract, we annouced the extraordinary success of quantum electrodynamics. How do we understand this 
success in the light of results obtained in this paper? In the following section, we discuss how a divergent series could 
possibly lead to a meaningful theory of quantum electrodynamics. 



V. A BRIEF MATHEMATICAL DIGRESSION: ASYMPTOTIC EXPANSION SERIES 



There is an imporatant class of series, known as asymptotic series, which frequently appear in physical problems. 
These series behave like a convergent series upto a certain number of terms but after that it behaves like a divergent 
series. This type of series is called asymptotic series [29 . |30j and is generally defined through a power series represe- 
tation of a function. The behaviour of an asymptotic series is very transparent in the following example (a version 
of Stirling's formula) given by Bender and Orszag (page 218 in [301] ) for the asymptotic series expansion of factorial, 
Nl = (Z - l) !, in powers of 1/Z . 



(Z - 1)! = {2^/Z) 1 ' 2 e- z Z z ( 1 + -Z- 1 + —Z- 2 - -^-Z^ - 571 Z' 4 + 
v ' y ' ' V 12 288 51840 2488320 

163879 z _ 5 5246819 z _ 6 534703531 z _ 7 4483131259 
209018880 + 75246796800 ~ 902961561600 ~ 86684309913600 + " 
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For example, for 0! , the terms get smaller for a while but the 15th term becomes larger than 0.01, and the 
35th is bigger than 10 10 . On the other hand, the 35th term for 9! is 10 10 /10 35 which is quite small but the 
175th term is bigger than one, and the 199th term is nearly 10 12 . This series has zero radius of convergence. 
A non-zero radius of convergence in 1/Z would also include some negative values around zero, and therefore, if 
the series for {Z — l)! converges for some large positive integer Z = {N + l) , it will also converge for some 
large negative integer Z = -(N + l). But(Z-l)!= ( - N - 2)\ is infinity. Therefore, the the 1/Z ex- 
pansion series of [Z— l)! for Z = N + l can not converge. The series is meaningful only as an asymptotic expansion. 

Mathematically, a function f(x) is said to have an asymptotic power series represetation if for all n, 

lim| /(3Q-EIW | = 

In other words, 

n 
i=0 

This means that the error in estimating the function is of the same order as the last term in the series. To explain, 
let us consider the following function, 



+ xt 

for real positive x and x — > 0. Since 



r°° e" 



1 2 2 (~Xt) k 

-l-xt + x 2 t 2 1 1 v ; 



1 + xt 1 + xt 

we have, 



N 

F(x) = Y / (-V k x k M + R N+ i(x) ; \Rn{x)\ = N\x n 

k=0 



The ratio of the two successive terms is 



x k k\ 

= xk 



a;(fc-i)(fc_ l)! 

This shows that the terms first decrease (since by assumption < x « 1) and then increase (when k > -). From 
this it follows that for a given value of x, there exists a best approximation. In other words, for a fixed value of x, only 
a definite accuracy can be achieved. However, the function defined by integral is well behaved and is non-analytic in 
x at x = 0. 

We have seen in previous sections that perturbation series of quantum eletrodynamics in coupling constant as well as 
in 1/Nf is divergent. Studies of the large order terms in perturbation series suggest that these series are probably 
asymptotic in nature. The fine structure constant of QED is 1/137 which is quite small, and therefore in the asymptotic 
series, one can consider terms up to a very large order. This can possibly explain the spectacular success of QED, in 
spite of fact that the series is divergent and asymptotic. 

In our view it does not make sense to look for any kind of singularity in coupling constant of a theory with asymptotic 
expansion and try to draw any big physical conclusion. There are attempts to give meaning to this kind of series 
through Borel summation but, again in our view, one does not achieve more than what one obtains through the 
summation of asymptotic series as explained above. 



VI. CONCLUSIONS 



In this paper, we have discussed several aspects of perturbation theory in quantum electrodynamics. Most impor- 
tantly, we have shown that the Landau singularity appears in the leading order terms in the 1/Nf expansion series. 
The restriction to leading order terms makes sense only when the series converges. We used Dysonian argument to 
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show that the convergence of the series leads to the instability of vacuum state in QED. This demonstrates that Lan- 
dau singularity reflects the instability of vacuum state. These problems can be avoided if the scries diverges. Divergent 
series as asymptotic series can provide physically meaningful results within some unavoidable errors depending on the 
value of the expansion parameter. The fine structure constant, which is the expansion parameter in QED, is small, 
and therefore, there is no surprise why QED is such a successful theory. 

Divergence of the perturbation series, suggested by Dysonian argument, is to save the vacuum state from the catas- 
trophic disintegration. But that is not the sole point of Dysonian argument. It also suggests non-analyticity of 
observables as a function of the coupling constant, e 2 when e 2 = ( similarly in 1/Nf for Nf = oo). Note that the 
non-analyticity in coupling constant seems to be in the infrared (IR) region. On the other hand, the non-analyticity 
of 1/Nf expansion series for Nf = oo is connected with Landau singularity which is in the ultrviolet (UV) region. 
It is not clear how the IR and UV regions are connected. However, we must remember that we are dealing with 
pathological situations where there are singularities and divergences. The main point of Dysonian argument, in our 
view, is that the non-perturbative QED is non-analytic and this behaviour can not be captured by a perturbative 
power series. More than fifty years of theoretical research has not been able to find a non-perturbative formulation 
which can remove the pathological aspects of QED. May be it is time for experimentalists to step in and look for ways 
to find non-analytical dependence of some physical observales on the coupling constant or flavours. 

VII. ACKNOWLEDGEMENTS 

I would like to thank here a large number of friends and colleagues with whom I discussed the problem of Landau 
singularity for almost a quarter century, and I can not possibly recall names of all of them. Special thanks are due to 
Prof. Spenta Wadia from whom I learnt about the Landau singularity and Dyson's proof of divegence of perturbation 
theory in the year 1983 when I was a graduate student in Tata Institute of Fundamental Research, Mumbai, India. 
I would also like to use this opportunity to thank Prof. N.V. Mitskievich, my M.Sc. thesis adviser and Prof. P.P. 
Divakaran, my Ph.D. thesis advisor, from whom I learnt various aspects of classical and quantum field theories. 



[1] L.D. Landau, in Niels Bohr and the Development of Physics, edited by W. Pauli (Pergamon Press, London, 1955); 

L.D. Landau and I.J. Pomeranchuk, Dokl. Akad. Nauk., SSSR 102 (1955) 489 
[2] D.V. Shirkov, Evolution of the Bogoluibov Renormalization Group, pp 25-58 in "Quantum Field Theory: A 20th Century 

Profile" Ed. by A.N. Mitra, Hindustan Book Agency and Indian National Science Academy, New Delhi 2000 
[3] J. Kogut, E. Dagotto and A. Kocic, Phys. Rev. Lett. 61, 2416(1988);Phys.Rev.Lett 60 (1988) 772 
[4] E. Dagotto, A. Kocic and J. Kogut, Phys. Lett. 232B (1989) 235 
[5] S. Kim, J. Kogut, and M-P. Lombardo, Phys. Lett. B502 (2001) 345 
[6] V. Miransky Nouvo Cimento 90A (1985) 149 

[7] M. Gockeler, R. Horsley, V. Linke, P. Rakow, G. Schierloltz, and H. Stiiben, Phys. Rev. Lett. 80 (1998) 4119 
[8] H. Gies and J. Jaeckel, Phys. Rev. Lett. 93 (2004) 110405 
[9] K. Langfeld, L. von Smekal, and H. Reinhardt, Phys. Lett. B362 (1995) 128 
[10] F.J. Dyson and A. Lenard, J. Math. Phys. 8 (1967) 423; 9 (1967) 698 

[11] EH. Lieb et al, The stability of matter: From atom to stars, Selecta of EH. Lieb, Ed. by W.E. Thirring, Springer- Verlag, 
1991 

[12] EH. Lieb et al, Phys. Rev. Lett. 79 (1997) 1785 

[13] J.G. Conlon et al, Commun. Math. Phys. 116 (1988) 417 

[14] E. Lieb and H.T. Yau, Commun. Math. Phys. 118 (1988) 177 

[15] J. Feldman, J. Magnan, V. Rivasseau, and E. Trubowitz, Euro-phys.Lett.24 (1993) 437; 24 (1993) 521 

[16] J. Froehlich, T. Chen and M. Seifert, "Renormalization Group Methods: Landau-Fermi Liquid and BCS Superconductors", 

|cond-mat/950 8063 
[17] R. Shankar, Rev.Mod.Phys. 66 (1994) 129 

[18] A. A. Abrikosov, L.P. Gorkov, and I.E. Dzyaloshinski, Methods of Quantum Field Theory in Statistical Physics, Prentice- 
Hall, INC., Englewood Cliffs, New Jersey, 1963 (Chapter-I, Sec.5, Page 40) 
[19] F.J. Dyson, Phys.Rev. 85 (1952) 631 

[20] Arkady I. Vainstein , "Decaying Systems and Divergence of Perturbation Theory", Novosibirsk Report, December 1964, 

reprinted in Russian with English translation by M.Shifman in proceedings of QCD2002/ArkadyFest. 
[21] Gerald V. Dunne, "Perturbative-Noperturbative Connections in Quantum Mechanics and Field Theory" , |hep-t h / 0207046 
[22] M. Azam, "Some Comments on the Divergence of Perturbation Series in Quantum Eletrodynamics" , Mod. Phys. Lett. 
A21 (2006) 1161 

[23] GM. de Divitiis, R. Frezzotti, M. Gagnelli, M. Masetti and R. Petronzio, The Bermions: An Approach to Dynamical 
Fermions from Negative Flavour Numbers, Nucl. Phys. B455 (1994) 274; |hep-lat/9507020| 



13 



[24] J. Rolf and U. Wolf, Running Coupling Constant for Wilson Bermions, Nucl. Phys. Suppl. 83 (2000) 899; |hep-lat/9907007| 
[25] S.J. Anthony, C.H.L. Smith and J.F. Wheater, Phys. Lett. 116B (1982) 287 

[26] I.M. Gelfand, R.A. Milnos, and Z.Ya. Shapiro, Representation of the Rotation and Lorentz Groups and Their 

Applications, Pergamon Press, 1963. 
[27] MA. Naimark, Linear Representation of the Lorentz Group, Pergamon Press, 1964. 
[28] S. Weinberg, The Quantum Theory of Fields, Cambridge University Press, 1995. 
[29] L. Sirovich, Techniques of Asymptotic Analysis, Springer- Verlag, New York-Heidelberg-Berlin, 1971 

[30] CM. Bender and S.A. Orszag, Advanced Mathematical Methods for Scientists and Engineers, McGraw-Hill Book Company, 
New York, 1978 



